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Although at temperatures T 3> Aqcd the quark-gluon plasma (QGP) is a gas of weakly interacting 
quasiparticles (modulo long-range magnetism), it is strongly interacting (sQGP) in the temperature 
range (1 — 3) Tc. One aspect of these interactions is the existence of many binary bound states of 
quasiparticles. Only qq ones have been so far directly seen on the lattice, for charmed and light 
quarks, but other attractive channels in qg,gg are likely to have them as well. It was argued in 
our previous paper that such bound states account for a significant part of the bulk properties 
such as density and pressure. Using the same model, we evaluate the energy loss dE/dx due to 
the ionization of these states. We found that it is substantial, but only in the narrow interval of 
temperatures T — (1.4 — 1.7)Tc. In contrast to that, we show that radiative and elastic losses are 
not likely to be modified much by binding, as the total density of color charges is close to what it 
is for weakly coupled quasiparticles. These distinctions would be important for understanding the 
energy dependence of jet quenching. 



Jet quenching is a sort of "tomography" of the 
prompt excited phase triggered in high energy heavy ion 
collisions. Even very hard jets are expected to lose some 
energy during their passage through the system, thereby 
providing information about the early stages of the col- 
lision. The quenching factor Q{pt) is defined as the ob- 
served number of jets normalized to the expected number 
of jets calculated in the parton model*. Experimentally, 
jet reconstruction in a heavy ion environment is very dif- 
ficult to achieve. Therefore, all currently reported results 
for jet quenching refer to the observed/expected ratio of 
the yields of single hadrons. Two-particle correlations 
have also been studied experimentally, confirming exis- 
tence of forward and backward (in azimuth) correlations 
expected from jets. 

In the early theoretical assessments of jet quenching [1] 
the mechanism considered was a jet re-scattering on 
quasiparticles in the Quark-Gluon Plasma (QGP). Tak- 
ing into account radiation effects leads to larger quench- 
ing [2] for hard jets. However, in dense enough matter one 
has to correct the radiative energy loss due to destruc- 
tive interferences (the so-called Landau-Pomeranchuck- 
Migdal or LPM effect), which modifies and reduces the 
effect, see e.g. [3] and references therein. For a recent 
brief summary see also [4] . 

One way to enhance the quenching in the prompt phase 
is through a synchrotron-like QCD radiation [5]. How- 
ever, this only takes place during a very short initial time 



*Effects due to initial state interaction (nuclear shadow- 
ing) are included. In other words, the parton distribution 
functions are nuclear rather than hadronic, and follow from 
lepton-nuclei experiments. Parton rescattering in nuclei at 
the origin of the so called Cronin effect is included in the 
expected yield. Final state interactions are excluded. 



in the heavy ion collision, when (and if) the gluons can 
be treated as a coherent classical field (the color glass 
condensate [6] or a set of QCD sphalerons [7] ) . 

Experimentally, a relatively modest quenching of 
rather high energy jets going through cold nuclear mat- 
ter has been first observed in deep inelastic scattering, 
see e.g. [8]. In contrast to that, the very first RHIC data 
at large pt have shown quite spectacular jet quenching, 
by about one order of magnitude. It implies that only 
jets originating from the surface of the nuclear overlap 
region reach the detector. Strong azimuthal anisotropy 
characterized by V2 = {cos{2(j))) is observed, close or even 
somewhat exceeding the "geometrical limit" [9] . 

Addressing the origin of this strong quenching, one 
may naturally start with a question: Is this quenching 
just proportional to the parton density involved, as so 
many mechanisms predict? The apparently near-absence 
of jet quenching at CERN SPS energies ''' with the fact 
that the relevant multiplicity dN/dy{y = 0) at these two 
energies is different by less than a factor of 2, it seems 
likely that a new mechanism of quenching opens up at 
RHIC energies. Future experiments will investigate this 
further, and at the time of this write-up, experimentalists 
are busy analyzing the latest ^/s = 62 GeV AuAu run at 
RHIC. 

In this letter we evaluate a contribution to jet quench- 
ing, induced by the "ionization" of the recently found 
binary bound states at T > T^. One motivation for that 
is that in ordinary matter, the QED dE/dx for charged 
projectiles with gamma factors in the range 7 =1-10^, is 
known to be dominated by the ionization losses. 



'One more argument, not emphasized enough, is the ratio 
of direct photons to tt decays: it is about 2-3 at RHIC while 
only 0.2 or less at the SPS. 
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static potentials [13]; (iii) quasiparticle masses [14]; and 
(iv) bulk thermodynamics [15]. For studies of binary 
states in A/" = A SUSY YM theory, in which a paramet- 
rically strong QGP-like phase is possible, see [17]. Last 
but not least, a liquid- like picture of matter also provides 
a natural explanation [10] for a successful description of 
collective phenomena at RHIC by ideal hydrodynamics. 

In this letter we will not go over the detailed argu- 
ments and calculations, which can be found in the above- 
mentioned papers. However and for definiteness we need 
to explain the parameters of a particular model used for 
the estimates to be carried below. Table 1 lists all binary 
attractive channels of quarks and gluons, indicating the 
number of states and their corresponding squared effec- 
tive dipole charge. 

For the effective (T-dependent) number of flavors, we 
will use 



Nf = 2 + e-"=/'^ 



= 120 MeV 



(1) 




fYfVi»t1*(*lYi I I I I I I I I I I I I I I I I I I I I I 



1 1.11.21.31.41.51.61.7 

FIG. 1. (a) The lines show twice the effective masses for 
quarks and gluons versus temperature T/Tc- Note that for 
T < 3Tc we have Mq > Mg. Circles and squares are the es- 
timated masses of the pion-like and rho-hke qq bound states, 
while the crosses stand for all colored states, (b) Density of 
various components of sQGP, normalized to T^, n/T^ versus 
the temperature T /Tc in units of the critical temperature. 
Circles and squares correspond to colorless mesons and col- 
ored bound states, respectively. The diamonds and crosses 
correspond to quark and gluon quasiparticles. The upper and 
lower lines are the total density of "states" , binaries and quasi- 
particles, and the total number of "charges". 



Although listed in the table, two channels which have 
smaller attraction than qq, namely qg^, qq, are ignored in 
estimates to follow. 

In Fig. la we show the (doubled) quasiparticle masses 
and those for bound states, as determined in a model used 
in [16] for global thermodynamics. For the purposes of 
this paper we need the particle densities, shown in Fig. lb. 
The highest T on that plot roughly corresponds to that 
at which qq, ggg states gets unbound. One can see that 
quark and gluon quasiparticles do not dominate the den- 
sity: the binary bound states do. It is however amusing 
to note, that the total number of "charges", defined as 
twice the density of composites plus those of quasipar- 
ticles, oscillate ^ around richarges/T^ « 4 in the whole 
interval. This is close to the total number for 16 mass- 
less gluons and 2A q + q {Nf = 2), the same within the 
uncertainties involved. 



channel 


rep. 


dipole factor Cd 


no. of states 


gg 
gg 


1 

8 


3 

9/4 


9. 

9. *16 


qg + qg 
qg + qg 


3 
6 


11/6 
4/3 


3c * 6s * 2 * Nf 
6c * 6, * 2 * Nf 


qq 


1 


4/3 


8, * Nj 


qq + qq 


3 


1 


4, * 3c * 2 * iV| 



TABLE L Binary attractive channels, the subscripts s,c,f 
mean spin, color and flavor, and Nf is the number of relevant 
flavors. 



Bound states in sQGP . Recently a radically new 
view of the QCD matter in the temperature range T =(1- 
3)Tc has emerged [10], in which the interaction between 
quasiparticles is strong enough to generate multiple bi- 
nary bound states of quasiparticles. This picture pro- 
vides a consistent description [11,16] of several previ- 
ously disconnected lattice observations, such as (i) bound 
states for charmonium and some light qq states [12]; (ii) 



*Note that a minimum around T = 1.2 Tc is where the ne- 
glected contributions of qq, qg% will peak. 



Thus, the parameters of our model for sQGP (fixed 
from the solutfon of the Klein-Gordon equation in [16] 
with lattice-based potentials and masses) posess an inter- 
esting but mysterious fine tuning, leading to a rather ac- 
curate duality between the total density of "charges" with 
a naive weakly coupled QGP made of massless quasipar- 
ticles. 

Photodissociation . We start by recalling some text- 
book results in QED. In the dipole approximation, the 
ionization cross section relates to the matrix clement of 
the momentum operator between the initial and final 
states, 



us [18], and subsequently by Bhanot and Pcskin [19] us- 
ing (4)to discuss charmonium disintegration by hadrons. 
A relativistic version of (4) with detailed balance built 
in was discussed by one of us [20] in the context of char- 
monium disintegration by gluo-effcct in the quark-gluon 
plasma. 

Our case is intermediate between (3-4): the bound 
states in the QGP have the effective Coulomb interac- 
tion at small r, which is however Debye screened at large 
r. One thus expect the threshold behavior of (4) and 
large energy behavior of (3). Let us also note that in 
the cases considered here, like for deuteron, no internal 
excitation possible as there is only one bound state. 
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(2) 



There are two cases in which the corresponding cross 
section has been analytically calculated: i. Coulomb 
systems, with a long-range interaction in the outgoing 
channel, and ii. systems bound by a short-range interac- 
tion, where the outgoing state can be approximated by a 
free plane wave. 

For Coulomb systems the cross section is known to be 
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(3) 



where the binding energy is A, the Bohr radius is a = 
1/mZa, the Gamow parameter is i^ = Za/v, with v 
being the outgoing velocity. The ionization energy is 
A = Z^a^'m/2. The remarkable feature of this cross sec- 
tion (due to long-range Coulomb final state interactions) 
is its finite value at the threshold lo — A, where the out- 
going velocity j/ — > oo with the last factor in (3) becoming 
finite 1/e^. In the opposite limit ui ^ A one finds i^ ^ 
and the cross section rapidly decreases (j{uj) « l/w^/^^ 

For short range interactions, an instructive example 
is the deuteron photodissociation. Using the deuteron 
wave function in a simple form '4'{r) « e~"''/r, where 
K^ — MA with M, A are the nucleon mass and binding 
energy, the photodissociation cross section is ^ 



a{u) 



Sna VA (w - A) 
~3 M 



3/2 



(4) 



In this case the cross section at threshold vanishes with 
the cube of the momentum following from (2), and peaks 
at w = 2A. 

In the QCD context the same effect has been discussed 
for the disintegration of heavy quarkonia by incoming 
gluons. This process was originally discussed by one of 



''We remind the reader that the effective charge in a dipole 
has in this case an additional factor of 1/2, since the neutron 
is not charged. 



Energy Loss for a high energy partons is estimated 
using the well known Weizsacker- Williams (WW) ap- 
proximation, which in QED views the Coulomb field of 
rapidly moving charges as a set of "equivalent photons" . 
Their number per frequency is 



, , „ a , f E\ duj 
n(Lu) = 2- In — — 

TT \ UJ J UJ 



(5) 



where the logarithm originates from the integral over the 
transverse momenta. In the QCD context this idea is 
known as "splitting functions" where one-kind of parton 
can split into two. At small relative energy the gluon-in- 
quark and gluon-in-gluon number per frequency is just 
(5) where a -^ ag^, with ^ — C(3) and S, = C(8) respec- 
tively. 

The ensuing energy loss per length dE/dz maybe writ- 
ten as 






The energy dependent mean-free path 

1 



A(c.) = 



p(7pb{u 



(6) 



(7) 



involves the parton-bound-state cross section apt and the 
density of bound states of mass M. At temperature T > 
Tc the latter is just {Ti — 1) 



p-M/T 

p = ^-^ = e-^/^ (MT/27r)3/2 

Arp 



(8) 



The integral in (6) leads dE/dx w \n{E/A) for relativis- 
tic particles, which is well known and well tested in QED 
in atomic ionization losses. 

The ionization cross section follows from the forward 
gluon-bound-state amplitude by the optical theorem, 
agb = ImMgb/w. For energies w « A, the forward 
amplitude is dipole dominated. Standard second order 
perturbation theory yields 

ImMgft == (g6|(d^* E^*) ttJ (h + A - lj) {d^^ E^^)\gb) (9) 

with the colored dipole moment for the bound state 



lAi 



{T^^T^)^^ 



(10) 



for two constituents each of mass m, respectively in the 
representation 1,2. The effective factor in cross section 
can be rewritten solely in terms of Casimirs 



1 



1 



Cd 



(Ci + C2) - -C_ 



(11) 



A comparison to our result written in a similarly approx- 
imate form yields the ratio 



l\E., 



2 E.g>' 1 

/S.Erad 9as Ugluons LM 



photo 



= (lo/L) 



(17) 



and it is listed in the Table 1. 

In the subsequent analysis we assume that the ionized 
outgoing constituents are free and non-relativistic with 
h « p^/tti. Averaging over the gluon spin and color 
yields 



{g\E^^E^^g)^^S^^S^^ 



(12) 



For Coulomb bound states with a deuteron-type wave 
function, straightforward algebra using (9-11) yield the 
ionization cross section 



CTgbiuj) 



52cd^(w-A)3/2 



(13) 



This result is in agreement with (4) modulo color factors. 
The temperature dependence of the inverse absorption 
length (A * Tc)^^ for a gluon of 4 different energies is 
shown in Fig. 2 (a). 

In terms of (13) the energy loss (6) is simply 



dE 
dz 
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(14) 



Using the density of binaries discussed above, and per- 
forming the integral over the energy of "equivalent glu- 
ons" we get our final results for the energy loss dE/dx 
for gluon jets as shown in Fig. 2b by the thick solid lines. 
(As usual, for quarks the results are 4/9 times smaller.) 
For comparison, we plotted (thin dashed lines) a con- 
tribution from elastic losses, using Bjorken's expression 



dE 
dx 



1 



C2 37ra^T^ln(3£;T/(2^^)) 



(15) 



where we used the Debye mass for ^ « 2T. Although 
this formula was derived for massless quasiparticles, due 
to the duality mentioned above it holds approximately 
for our model with binaries as well. We note that al- 
though the magnitude is comparable, the T-dependence 
of dE/dx due to these two mechanisms is quite different. 

Radiative losses . A comparison to the radiative 
losses can be better done not in terms of dE/dx but in 
total losses after a path of length L, since the LPM effect 
makes dE/dx ^ L [3]. Using the same simplifications 
as discussed in [4], we can write radiative energy loss on 
gluons as 



AE, 



2lTTa\ 



rad 



ngiuonsL^ ln(£'/^) 



(16) 
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FIG. 2. (a) The inverse absorption length for a gluon in 
sQGP, in units of Tc, versus the temperature T/Tc in units of 
the critical one. Circles, squares, diamonds and crosses cor- 
respond to a gluon energy oj = 0.5, 1, 1.5, 2 GeV, respectively, 
(b) Gluon energy loss dE/dx in GeV/fm versus the tempera- 
ture T /Tc (in units of the critical one) . The thick solid lines 
are for the "ionization" losses, while the thin dashed lines are 
for the elastic losses. In each set the three curves from top to 
bottom are for a gluon with energy 15, 10 and 5 GeV. 

with /o — (.5 — Ifm). (The sum above is over the 
binary species and M is the quasiparticle mass.) The 



radiative losses are thus larger than ionization ones if the 
length the jet passes in matter is larger than ^o- 

However, one can easily see that this difference is 
hardly important in practice. Indeed, already from the 
magnitude of the experimental jet quenching Q{pt) ~ 
0.1-0.2 and the rather steep slope of the pt spectra in 
the relevant pt region, it follow that a loss of about 
AE/E « 1/5 is sufficient at RHIC conditions to get 
a suppression of the order of that observed, making 
such jets invisible anyway. As one can see from Fig. 2, 
such losses will take place for a gluon (or quark) pass- 
ing through only L « 0.5 fm (or 1 fm) of matter with 
T « l.QTc « 300 MeV. For such length both quenching 
mechanisms are about equally important. 



Discussion. It is well known that partonic jets of 
quarks and gluons, promptly produced in RHIC colli- 
sions, can be quenched via gluonic radiation. We have 
shown that additional "photodissociation" mechanism of 
quenching can appear, provided there are binary bound 
states at T > Tc. The energy loss is few times A the 
binding energy of the gg, qq or gq per collision. We eval- 
uated the total energy loss due to it, and found that it is 
comparable to the radiation loss in magnitude for lengths 
L w 1 fm. 

Due to the "duality" between the total density of 
charges in our model and the naive gas of massless quarks 
and gluon noticed above, the radiative energy loss is the 
same for both (strong) sQGP and naive (weak) wQGP 
models. The "photodissociation" we have discussed, on 
the other hand, is present only in the sQGP scenario, 
and is only there in a restricted interval of temperatures 
Tc < T < 1.7 Tc w 300 Mey. Such T are reached at 
RHIC but not at the SPS. This distinction, as well as 
very restricted interval of T where new mechanism ex- 
ists, will help us to decide which scenario is the case. 

Another important distinction between the two mech- 
anisms of quenching is that radiation remains in the for- 
ward cone, and so it can be found there experimentally. 
The energy lost by ionization of binaries obviously is sim- 
ply dissipated in the heat bath. 
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